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We report the results of an inelastic X-ray scattering study of the lattice dynamics in the para¬ 
electric phase of the antiferroelectric lead hafnate PbHfOa. The study reveals an avoided crossing 
between the transverse acoustic and transverse optic phonon modes propagating along the [1 1 0] 
direction with [1 -1 0] polarization. The static susceptibility with respect to the generally incom¬ 
mensurate modulations is shown to increase on cooling for the entire T-M direction. We consider 
different approaches to the data analysis that correspond to different models for the temperature 
evolution of the dynamic susceptibility function. A number of similarities and differences between 
the lattice dynamics of PbHfOa and PbZrOs are described. 

PACS numbers: 77.22.-d, 77.65.-j, 77.90.-|-k 


I. INTRODUCTION 

Perovskite ferroelectrics and similar crystals find nu¬ 
merous applications in technology^. These materials 
form a rich test ground for assessing the phenomena 
associated with the structural phase transitions. Since 
the introduction of the soft mode concept by Cochran^, 
and Anderson^ in the context of the ferroelectric transi¬ 
tion in BaTiOs, the topic of the lattice dynamics in per- 
ovskites attracted considerable attention. A number of 
thorough reviews were published before 1990^“— that con¬ 
sidered the results of the soft mode spectroscopy in var¬ 
ious perovskites. More recently an inelastic X-ray scat¬ 
tering (IXS) study of the critical dynamics in PbZrOa- 
like antiferroelectrics revealed that the critical dynamics 
in the paraelectric phase of PbZrOa is compatible with 
the one expected for the crystal with an incommensurate 
phase transition^. The key feature is a flat transverse 
acoustic (TA) branch, which shows an almost homoge¬ 
neous softening along a large part of the P-M (wavevector 
q = (9i 0)) directional^. 

In this paper we report the IXS study of a similar 
crystal, lead hafnate, PbHfOgi^. It is different from 
PbZrOs by the presence of an intermediate phase be¬ 
tween the cubic phase and low-temperature antiferroelec¬ 
tric phasei2r2i. The structure of the intermediate phase 
(or phases, according to some reportsi^) is currently un¬ 
der debate. The current study covers the lattice dynam¬ 
ics in the high-temperature cubic phase. 

In section El we specify the experimental methods, in¬ 
cluding sample preparation and IXS measurements. In 
section Ell we specify the basic model used to fit the IXS 
spectra. The experimental results are covered in sections 
m and El dedicated respectively to the analysis of dis¬ 
persion curves at high temperatures and to the analysis 


of the temperature evolution of the signal. A compact 
summary of the results and the data treatment is given 
in section ED The last two sections are dedicated to the 
analysis of possible interpretations and for the conclu¬ 
sion. 


II. EXPERIMENTAL METHODS 


A. Sample 

The single crystals of PbHfOa were synthesized by 
means of spontaneous crystallization from the high tem¬ 
perature solution in a Pb304-B203 solvent. The compo¬ 
sition of the melt used in our experiments was as follows: 
2.4 mol% of PbHf03, 77 mol% of PbO (re-counted to 
Pb304) and 20.6 mol% of B2O3. The Pb304 was used 
instead of PbO to avoid coloration of the as-grown crys¬ 
tals caused by oxygen deficiency. The crystallization was 
carried out in a platinum crucible covered with a plat¬ 
inum lid under conditions of a low temperature gradient. 
After soaking at 1473 K for 24 hours the melt was cooled 
down to 1200 K at a rate of 3.5 K/h and after decantation 
the furnace was cooled to room temperature at a rate of 
10 K/h. As-grown PbHfOs single crystals were etched 
in diluted acetic acid to remove residues of the solidified 
flux. The sample for the IXS measurements was a 40 by 
40 by 600 micron size needle. In accord with the litera¬ 
ture, the cell parameter of the cubic phase was found to 
be a = 4.14 at T = 473 K. 
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B. IXS measurements 

The experiment was performed at ID28 beamline of the 
ESRF (European Synchrotron Radiation Facility), using 
the Si (9,9,9) setup at E = 17.794 keV which provides an 
overall energy resolution of about 3 meV— . We studied 
the phonons propagating along the [1 0 0] and [1 1 0] 
direction with polarization (direction of the atomic dis¬ 
placements) lying in the plane defined by these two vec¬ 
tors. Thanks to the fact that the spectrometer hosts 9 in¬ 
dependent analyzer-detector pairs, the spectra at points 
away from these directions but within the correspond¬ 
ing plane were also recorded. The sample was glued to 
quartz capillary. The temperature was controlled by the 
ESRF mini heat blower from 473 K (just at the bottom 
of the stability region of the cubic phase) to 773 K with 
an accuracy of about 2 K. 

III. IXS SPECTRA MODELING 

The modeling of the spectra was done following the 
standard approach for which the one-phonon scattering 
function for a phonon with wavevector q in branch j is^^ 

Sj{Q,u}) = |Gj(Q,q)|2Fj [uj,ujj{q),T ], (1) 

where Gj(Q,q) - the inelastic structure factor and 
Fj[uj,ujj(q),T] - the response function, hio - the energy 
transfer in the one-phonon scattering event and Wj(q) 
describes the dispersion curve. The magnitudes of the 
inelastic structure factorsi^ were treated in our analy¬ 
sis as variable parameters. The model response function 
corresponds to the dynamic susceptibility of a damped 
harmonic oscillator, given byi^ 

xSho = - iujT)-\ (2) 

where T is a damping constant. The resulting expression 
for the response function isi^ 

Fj[uj,ujjiq),T] = 

^ _ T) _ 

1 - exp{-hjj/kT) [uj^ - a;|(q, r)]2 -|- uj^r‘^{q,T)' 

The IXS signal was modeled by a numerical convolu¬ 
tion of the sum of the scattering functions ([T]) (DHO) 
with the resolution function of the IXS spectrometer. 
At a number of occasions we also added the resolution- 
limited central peak (CP). The resolution function was 
modeled by a pseudo-voigt lineshape with independently 
determined parameters. 

IV. RESULTS: PHONON DISPERSION AT 
HIGH TEMPERATURE 

We start from the general overview of the dispersion 
curves at high temperatures. At T=773 K the two 
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FIG. 1. Phonon dispersion curves of PbHfOs at T=773 K. 
The solid lines show a fit by the rigid ion model. 


phonon resonances are well resolved. This allows us 
postponing the discussion of the less straightforward low 
temperature data treatment process to the next sections. 
The dispersion curves are shown in Fig. [T] 

The key feature is the TA (transverse acoustic) branch, 
which displays a relatively weak g-dispersion in T-M di¬ 
rection (q = {q,q,0)) and appears relatively flat for 
q > 0.075 (hereafter q is expressed in reciprocal lattice 
units, r.l.u.). This is similar to the result of Refs. Handle 
for PbZrOa. Approximately at q = qAC ~ 0.075 there is 
a rapid change in the group velocity. The linear slope of 
the transverse optic (TO) branch at q > qAC appears to 
be approximately the same as the linear slope of the TA 
branch for q < qAc (which corresponds to the transverse 
speed of sound). So, they can be viewed as continua¬ 
tions of each other. As it appears from Fig. [T] there 
is an avoided crossing (also referred to as anti-crossing) 
between the TA and TO branches. In perovskites this 
phenomenon is not uncommon. In particular the stud¬ 
ies of KNbOa^ and BaTiO^i^ show a similar behavior. 
The shell model fit in KNbOa suggests^^ that the flat 
part of the TA branch largely corresponds to the mo¬ 
tion of the ferroelectric-active niobium ion. By analogy 
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one may suggest that the motion in the flat part of the 
TA branch corresponds to the motion of the ferroelectric- 
active ion. A similar assumption was made in for¬ 

mulating the simple model (“polarizability model”) for 
PbZrOa and PbHfOa in Refj^^. 

We performed the analysis of the motion in the sim¬ 
plest possible harmonic model - the rigid ion model 
(RIM). For the perovskite structure this model is de¬ 
scribed in Refs, and[2l|. The solid lines in Fig. [T] cor¬ 
responds to the least squares fit by RIM with the follow¬ 
ing parameters: longitudinal force constants Ai=20.15, 
A2=289.1, A3=-11.22, transverse force constants Bi=- 
3.944, B2=-63.09, B3=3.595 (in units e^/2ua where e - 
electron charge and Uo=(4.15 A)^ - unit cell volume), 
charges of lead and hafnium Z\ = 1.7e and = 4.07e. 
The fit by RIM indicates that the TA branch in the F-M 
direction, for q > qAC-, corresponds mainly to the motion 
of lead and hafnium with displacements (not normalized 
by the square root of the mass) of lead being about 2 
times larger than the ones of hafnium. This supports 
the assumption of mainly lead motions for the flat TA 
branch. In principle, a more reliable analysis could be 
possible with more detailed harmonic models, such as 
the shell model, which accounts for the electronic polar¬ 
izability of the ions. This type of analysis is out of scope 
of the current paper, in part due to the absence of reliable 
data on the high-energy modes in the paraelectric phase 
of lead hafnate, which are required for such analysis^S. 


The static one-phonon susceptibility is defined by the 
matrijii 


ySta‘(q) = /3(Qj(q)Qy(-q)), (6) 

where /3 = {kT)~^, () - ensemble average and Qj(q) - 
normal mode coordinates. In the case of susceptibility 
with respect to the ionic displacements within a particu¬ 
lar phonon mode, the diagonal terms in Eq. m are given 
(in the high-temperature approximation) byi 


Af"‘(q) = 


1 


( 7 ) 


On the other hand, this static susceptibility is propor¬ 
tional to the integral of the scattering function m over 
all energy transfers 


^Stat 

^3 


(q) = 


/3 


|GAQ,q)l 


5r(Q). 


2 3 


( 8 ) 


Finally, as it might be the case of PbHfOs, the crystal 
may have degrees of freedom that are not of harmonic 
nature. In particular, they can be slowly-relaxing quasi¬ 
spin ionic shifts. In this case the dynamic susceptibil¬ 
ity will not be of the form ([2]), but the relations o,® 
and ® will be still relevant, provided that Qj(q) and 
Gj(Q, q) refer to the appropriately defined normal mode 
coordinates and structure factors. 


V. RESULTS: TEMPERATURE EVOLUTION 
OF THE IXS SIGNAL 


A. Relations for susceptibilities 


B. Dielectric permittivity and the approach for 
analyzing the temperature dependences of finite-g 
susceptibilities 


The inelastic and quasielastic scattering in the high- 
symmetry phase is used in this study to gather the infor¬ 
mation on the susceptibility of the crystal with respect to 
various sorts of ionic displacements that break the high 
symmetry^. It is often the case that on approaching the 
structural phase transition of a particular type, the static 
susceptibility with respect to the corresponding patterns 
of ionic displacements grows and even diverges at the 
transition point. The simplest case, which we assume 
throughout the paper, is when the transition is second 
order like and the static susceptibility obeys the Curie- 
Weiss (CW) law 


(4) 

where C is the CW constant. The static susceptibility 
in this formula should not be confused with the dynamic 
susceptibility. The connection between these two suscep¬ 
tibilities is given by the Kramers-Kronig relation^ 


= Re [xDy"(a; = 0)] = - H 

TT 7-00 W' 


( 5 ) 


Due to the finite resolution effects, it is not possible 
to assess small-energy phonons at the Brillouin zone cen¬ 
ter by the IXS technique. Published studies of the tem¬ 
perature dependence of the dielectric permittivity sug¬ 
gested that Xoieiectric obeys CW law. The CW tem¬ 
perature was reported as yShirane _ ^23 K by Shirane 
et ah^ and as = 378 K by Samarai^. In 

these two studies the CW law was determined by ana¬ 
lyzing the dependence of the dielectric permittivity e(T) 
in the range T). < T < Tc + AT, where AT is about 90 
K. The recent report by Bussmann-Holder et al^ (45 
years after Samara’s work), where e(T) was measured in 
a broader temperature range, shows that the situation 
can be more complex. Although it appears that their 
data for T^ < T < Tc + AT are roughly consistent with 
To = jjSamara _ gyg (data from Ref. [ 2 ^ would have 
given To = 404K assuming the same Tc as in Samaras pa¬ 
per), for higher temperatures, T > Tc + A, the e(T) obeys 
CW law with Tq = To^^^^ann ^ ^ -g 

reported as 451 K, but the private communication with 
the authors indicated there was a misprint). It is worth 
mentioning that = 441 K is close to the tem¬ 

perature of the phase transition to the low-temperature 
antiferroelectric phase [Tai-A 2 in Ref. 1^ . 
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In order to compare the temperature dependences of 
different kinds of susceptibility we define a function de¬ 
scribing the relative rate of growth as 


R{T) = 


1 

xStat(2^) dr 


(9) 


In terms of this function we may state that the relative 
rate of growth of e(T), for Tc < T <Tc + AT, is numeri¬ 
cally similar to the one for described by CW law 

with To « The function R{T) can be useful in 

comparing the temperature dependences of susceptibili¬ 
ties when the absolute scale of is unknown, as it 

is the case in scattering experiments. In the case that the 
susceptibility can be described by x®*®'* = C/{T — Tq)", 
the function R{T) = nj(T — Tq) contains the information 
on both To and power n. 


C. Modeling the temperature dependence of the 
signal at q = (0.1, 0.1, 0) 

Fig. [5]shows the temperature dependence of the signal 
at q = (0.1,0.1,0). The two panels (left and right) show 
the same data fitted by two different models that will be 
discussed in this section. It is evident that on temper¬ 
ature decrease the scattering at small energy transfers 
gets stronger. The overall shape of the low-energy part 
of the spectrum (everything excluding the second, high- 
energy phonon peak) gradually transforms from the un¬ 
derdamped, two-peak DHO lineshape to the single peak 
centered at zero energy transfer. At the lowest temper¬ 
ature the peak appears to be almost resolution limited. 
Neither the energy nor the damping constant of the sec¬ 
ond (TO) phonon appears to be temperature dependent, 
as deduced from visual inspection of the spectra and from 
the fits obtained by different methods that we employed. 
The data allows different approaches for their treatment. 
In this section we discuss 3 scenarios, in which we ac¬ 
count for the temperature dependence of the signal by 
(1) temperature-dependent TA branch, (2) temperature- 
dependent central peak and (3) a combination of the two. 


1. Scenario 1. Temperature-dependent (soft) TA branch. 

This scenario is the formally simplest model function, 
consisting of two DHO contributions corresponding to 
the TA and TO phonons respectively. We set up the 
fitting process so that all the spectra, corresponding to 
a particular q (currently: q = O.I) and the whole set of 
temperatures, are fitted together. This allows improving 
the reliability of the optimization problem by setting up 
constraints on the model parameters that are not possible 
with independent fitting of different spectra. We set up 
the following constraints: 

• The energy and the damping constant of the TO 
phonon is temperature-independent. 


Scenario 1, q=0.1 Scenario 2, q=0.1 

T=773 K, q=0.1 T=773 K, q=0.1 



T=693 K, q=0.1 T=693 K, q=0.1 




T=603 K, q=0.1 T=603 K, q=0.1 



T=473 K, q=0.1 T=473 K, q=0.1 




Energy (meV) Energy (meV) 


FIG. 2. Temperature evolution of the IXS spectra at Q = 
(2.1, 0.1,0). The left and right panels show the same data 
fitted using the two different models described in the text. 
The data are normalized by the value |GtoPT where |GtoP 
is the squared absolute value of the TO phonon structure 
factor, as obtained by the fit, and T - temperature. 
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• The ratio between the structure factors of TA and 
TO phonons is temperature-independent. 

The first constraint is motivated by the fact that we 
did not observe any notable difference in the relevant val¬ 
ues when fitting the different spectra independently. The 
second constraint is due to the assumption that the struc¬ 
ture factors do not change with temperature. The rea¬ 
son for this is that the pattern of the ionic displacements 
does not change with temperature and the temperature 
evolution manifests itself only in the temperature depen¬ 
dences of energy and damping constant of the phonon. 
In principle, the second constraint could be reformulated 
by fixing the absolute values for both structure factors. 
However, in this experiment we were unable to secure 
the same intensity scale for scans at different tempera¬ 
tures and different momentum transfers (this is generally 
challenging), so we were able to fix only their ratio. The 
result of the described fit is shown in the left panel of 
Fig.l 

The temperature dependence of the TA phonon energy 
and of the damping constant are shown in Fig. [31 The 
energy gradually decreases while the damping gradually 
increases. The function Eq. © stops displaying the 
two-peak shape when F^ > u}j'/2, which corresponds to 
overdamping of the phonon. From the dependences in 
Fig. |3|it thus follows, that almost in the whole studied 
temperature region the response is underdamped and the 
apparent merging of the TA Stokes and anti-Stokes peaks 
already at T = 603 K is a resolution effect. Overdamping 
is not excluded at the lowest studied temperature. 

We obtain the model expression for the frequency of 
the TA phonon assuming that the corresponding suscep¬ 
tibility follows the CW law by combining Eqs. © and 
©: 


^ta = ^a(T-rTA), (10) 

where a is a constant. Here we presume that the TA 
branch corresponds to the motion of the order parameter 
relevant to the phase transition characterized by modu¬ 
lation with a particular wavevector (in the current case 
q = 0.1). 

The corresponding least-squares fit is shown in Fig. |3| 
as a solid line. This results in value = 295 ± 115AT. 
The corresponding relative rate of growth R{T) (Eq. ([3])) 
is thus comparable to the one of e(T). 


2. Scenario 2. Temperature-dependent central peak. 

In this scenario, in addition to the TA and TO phonon 
resonances, a resolution-limited central peak is present. 
The constraints are the following: 

• The energies and the damping constants of the TA 
and TO phonons are temperature-independent. 


Scenario 1, q=0.1 



Temperature (K) 


FIG. 3. Temperature dependence of the TA phonon energy 
and damping constant for Q = (2.1, 0.1,0), obtained by the 
treatment in accord with the Scenario 1 described in the text. 
Error bars correspond to the 95% conhdence bounds. The 
Curie-Weiss fit to the phonon energy dependence is shown by 
solid line and the corresponding estimate of the Curie-Weiss 
temperature To is shown by a dashed line. 


• The ratio between the structure factors of TA and 
TO phonons is temperature-independent (as in Sce¬ 
nario 1). 


The fit by this model is shown in the right panel of 
Fig. [3] The simplest interpretation of the central peak is 
when we consider it as being due to an independent de¬ 
gree of freedom. The response of this degree of freedom 
to the external force is of relaxational character with a 
relaxation time (which determines its energy width) suf¬ 
ficiently longer than the period of the lattice vibrations. 
This relaxational degree of freedom may correspond to 
the hopping of ions between the wells of the multi-well 
local potential. We presume that this relaxational de¬ 
gree of freedom may be characterized by a susceptibility 
described by the CW law. Combining equations O and 
dSj) the temperature dependence of the integrated central 
peak intensity can be written as 




|g.-(Q,q)P g 

/3 T- TgCP ■ 


( 11 ) 


Fig. [4] shows the temperature dependence of the cen¬ 
tral peak intensity normalized to the absolute value of 
the square of the TO phonon structure factor, which is 
assumed to be temperature independent. The solid line 














6 


Scenario 2, q=0.1 



Temperature (K) 


FIG. 4. Temperature dependence of the central peak intensity 
for Q = (2.1, 0.1, 0), obtained by the treatment in accord with 
the Scenario 2 described in the text. Error bars correspond 
to the 95% confidence bounds. The Curie-Weiss fit is shown 
by solid line (fitted value of CW tempearture is To = 436 ± 36 
K. Dashed line shows CW fit with Tq fixed at value 378 K. 
The inset shows the temperature dependeces for the relative 
rates of growth corresponding to the curves in the main plot. 


shows the least-squares fit by Eq. CD- The fit gives a 
value of = 436±36 K. This value is numerically simi¬ 
lar to jjBussmann^ jg substantially larger than Tg 
The data was also fitted by a CW law with Tq fixed at 
jjSamara. result is shown by the dashed line. An ad¬ 
ditional fit by a formula equivalent to (ED, but with the 
denominator (T — replaced by (T — where 

n is an additional free parameter, gives Tq = 285 K and 
n = 2.4. This fit (dotted line) appears slightly more con¬ 
sistent, but also not completely adequate. None of these 
fits is really satisfactory. 

The inset of Fig. |4] shows the relative rates of growth 
(defined by Eq. ([9])), by the same line types, as in the 
main plot. Near the transition temperature (about 470 
K) R{T) for the susceptibility associated with the central 
peak (either dotted or solid lines) is larger than the one 
for e(r) (dashed line). 


3. Scenario 3. Both the TA phonon and the eentral peak 
are temperature-dependent. 

This scenario is a hybrid of Scenario 1 and Scenario 
2. The temperature dependence is admitted in the TA 


phonon energy and damping and in the magnitude of the 
central peak. The constraints are the same as in the Sce¬ 
nario 1. The results of the fit are shown in Fig. [5^. The 
temperature dependence of the TA energy and damping 
are weak (Fig. [5lo): the Curie-Weiss fit to the tempera¬ 
ture dependence of the energy gives Tq < 0 K. The cen¬ 
tral peak intensity is treated in the same manner as in 
Scenario 2 (Fig. [^) resulting in Tq = 447±9.5 K. Appar¬ 
ently, the less constrained Scenario 3 gives qualitatively 
the same results as the more constrained Scenario 2. The 
temperature dependence of the intensity is more consis¬ 
tent with the CW law than within Scenario 2. On the 
other hand, the relative rate of growth is, again, larger 
than the one for e(T). 


4- Comparison of the fit statistics for q—0.1 

From the visual inspection of Fig. [2] it is evident that 
Scenarios 1 and 2 produce equivalently convincing fits. 
Apart from the visual inspection, a common way of eval¬ 
uating the “fit quality” is by comparing the value of the 
root mean square error for different models, defined as^^ 


R = 


(/Exp_jModeb2 

^i=l 57 


N -M 


( 12 ) 


where - the data points, _ corresponding 

model calculations, at - expected variances for data 
points (size of the error bars), N - number of data points, 
M - number of variable parameters. It is the same pa¬ 
rameter which is minimized during the fit process. In Sce¬ 
nario I the value of i?si = 1-4, in Scenario 2 i?s 2 = 1.44. 
Consequently, this estimate of the “fit quality” assigns 
only a minor favor to Scenario 1 over Scenario 2. Sce¬ 
nario 3 has more independent variables than Scenarios 
1 and 2; in fact the models of Scenarios 1 and 2 are 
the limiting cases of the model in Scenario 3. The aim 
of Scenario 3 can therefore be, at the best, to examine 
what the data fitting algorithm will “prefer” when pro¬ 
vided with increased flexibility. Apparently it “prefers” 
Scenario 2. However, one shall expect it to “prefer” as 
well Scenario 1 when the starting point in the data fit¬ 
ting algorithm would be chosen accordingly, because the 
corresponding local minimum in the objective function is 
deeper for Scenario 1 {Rsi < Rs 2 )- 

The conclusion of this section is that the temperature 
dependence of the IXS spectra at q = 0.1 can be de¬ 
scribed equally well by the two different models (Scenar¬ 
ios 1 and 2) that allow a reasonable conditionality of the 
optimization problem. The “physical relevance” of these 
two models shall be assessed at a later stage, in con¬ 
nection with the results obtained at different wavevector 
transfers. 
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Scenario 3, q=0.1 

a) b) 



FIG. 5. Fit results within Scenario 3 for data at Q = 
(2.1, 0.1,0). (a) Temperature evolution of the spectra and 

the corresponding fits. Error bars correspond to 95% confi¬ 
dence bounds, they are smaller than symbols, (b) Tempera¬ 
ture dependence of the TA phonon energy (circles) and damp¬ 
ing constant (triangles). The Curie-Weiss fit of the energy is 
represented by the solid line, (c) Temperature dependence of 
the central peak intensity and the corresponding CW fits with 
variable To (solid line) and To = = 373 K (dashed 

line). 


D. Modeling the temperature dependence of the 
signal in the remaining part of F-M direction 

When dealing with the data recorded at Q = 
(2.1, 0.1,0) we tried the three scenarios. Although Sce¬ 
narios 1 and 2 have different number of independent fit 
parameters, they both cannot be further reduced by im- 


q _0.1 0.2 0.3 0.4 0.5 

Scenario 1 1.40 1.23 1.34 1.29 1.27 
Scenario 2 1.44 1.38 1.64 1.48 1.36 

TABLE I. The root mean square errors characterizing the 
fits of the data along the F-M direction using the two fitting 
scenarios, as described in the text. 


posing more constraints without losing the ability to give 
a “good fit”. Scenario 3, on the other hand, is over¬ 
parametrized which makes, in general, the fitting prob¬ 
lem less well conditioned. Nevertheless, it provided rea¬ 
sonable results for the IXS data at 5 = 0.1. The spectra 
at larger q are in general less informative, meaning that 
the peaks are less well resolved. This makes Scenario 
3 practically inapplicable - the fitting problem becomes 
ill-conditioned and gives irregular results. For the other 
g-points along the T-M direction we therefore used only 
the first two scenarios. The results of the fits are shown 
in Appendix lAl 

Generally, at larger q the TO phonon resonance is less 
resolved than at q — 0.1. For g = 0.3 and larger wavevec- 
tors it becomes overdamped. The visual appearance of 
the fits, as in the case of g = 0.1 does not allow favoring 
one scenario over the other. The result of the numeri¬ 
cal evaluation using dUD is shown in Table U for all the 
points along the F-M direction. In all cases the root 
mean square error is slightly less for Scenario 1. We do 
not claim, however, that this difference can be used to 
decide which scenario is more likely to be correct. 

We attempted to bring more insight into the problem of 
distinguishing the two scenarios by considering the tem¬ 
perature evolution of the spectra at wavevectors, which 
are slighlyt away from the F-M direction. In particu¬ 
lar, we considered the point q = (0.12,0.21,0), which 
is, on one hand, close enough to the F-M direction to 
be influenced by the same physical effects and, on the 
other hand, is far enough to have larger energy of the TA 
phonon to provide a better distinction of its temperature 
dependence. This analysis was, however, inconclusive. 

The temperature dependence of the TA phonon dis¬ 
persion curve along the F-M direction, as determined ac¬ 
cording to the Scenario 1, is shown in Fig. [51 The energy 
decreases monotonically for g > 0.1 at all temperatures. 
At large g, the fit with Scenario 1 gives slightly modi¬ 
fied energy for the phonons at T=773 K, as compared to 
the energy determined by independent fitting of just the 
high-temperature data. The independently determined 
high-temperature dispersion curve from Fig. [T] is shown 
in Fig. iniby a solid line. On the other hand, at g < 0.2 
Scenario 1 is consistent with the result of the indepen¬ 
dent fit. For Scenario 2 the situation is different: at 
g > 0.3 the solid curve (independent high-temperature 
fit) is within the error bars of the points obtained within 
Scenario 2, while for g = 0.2 the curve is substantially 
higher than the point from Scenario 2. 
































FIG. 6. Temperature dependence of the TA phonon energy 
along the F-M direction, as obtained by fits according to Sce¬ 
nario 1. Solid lines connecting symbols are guides to the eye. 
Smooth symbol-free solid line - calculation by the rigid ion 
model from Fig. [T] The dashed line shows the TA phonon 
dispersion, as determined according to the Scenario 2. 


E. R-point 


The spectra for the R-point (Q = (3.5, 2.5,0.5)) are 
shown in Fig. [7^. At this point we took the spectra 
at the same temperatures as the spectra for the F-M 
direction, except the lowest temperature (473 K). The 
spectra contain the Stokes and anti-Stokes resonances of 
one phonon and the central peak. 

The independent fits for these spectra admit a weak 
softening of the phonon resonance - about 1 meV on 
cooling from 773 K to 513 K. The energy and damping 
constant of the phonon in this region are estimated to be 
9 and 10 meV, respectively. The admitted level of soften¬ 
ing, relative to the average energy (9 meV) is compara¬ 
ble to the one suggested by the analysis within Scenario 
3 for the TA phonon at q = (0.1,0.1,0) (see Fig. [^d). 
Neglecting this possible weak softening, we performed 
the fit of the spectra at R-point similarly to the fits in 
Scenario 2 for the F-M direction, assuming the energy 
and damping constant of the phonon to be temperature- 
independent, the result is in Figl^. The temperature 
dependence of the central peak intensity (normalized by 
the squared modulus of the phonon structure factor) is 
shown in Fig. [7)3. We provide two CW fits of the inten¬ 
sity. The first fit (solid line), with variable CW tempera¬ 
ture, yields Tq = 500 ± 25 K . This temperature is above 
the transition temperature. For comparison we provide 
a fit with Tq fixed at the transition temperature (dashed 
line). The temperature dependence of the central peak 
intensity at the i?-point is inconsistent with the CW law 
with critical temperature below the transition tempera¬ 
ture. This renders the assumption that the central peak 
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1=773 K, R-point T=693 K, R-point 





Temperature (K) 

FIG. 7. (a) Temperature dependence of the IXS spectra, ob¬ 
tained at the R-point (Q = (3.5, 2.5,0.5)), and the results 
of the corresponding fit. (b) Temperature dependence of the 
central peak intensity and the Curie-Weiss fit with variable 
To (solid line) and with To fixed at 470 K (transition temper¬ 
ature) (dashed line). With variable CW temperature the fit 
gives To = 500 K, which is larger than the temperature of the 
transition. 


at the R-point can be directly associated with a suscep¬ 
tibility following CW law inconsistent. A more complex 
model is required for the consistent description of the 
data. 


VI. SUMMARY OF OBSERVATIONS 

In the previous section we gave thoroughly illustrated 
results of the IXS measurements and of their basic analy¬ 
sis. This analysis was limited to (a) fitting the IXS spec¬ 
tra by model functions, (b) fitting the dispersion curves 
by the rigid ion model and (c) fitting the temperature 
dependences by the Curie-Weiss law in an appropriate 
form. Here we summarize the findings. 
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1. F-M direction, transverse displacements 

1. At high temperatures, there is a relatively flat TA 
branch, with a sharp change in the phonon group 
velocity at q = qac- At this point there is an 
avoided-crossing between TA and TO branches. 

2. The TO branch, at q > qAC, is temperature- 
independent, at least within the accuracy of our 
data. 

3. On cooling down, there is an intensity increase of 
small-energy-transfer-scattering (SETS). By SETS 
we mean the scattering due to everything except 
the TO phonons. 

4. The data admits different models for the SETS 
growth. 

5. In Scenario 1 the SETS growth is due to the TA 
branch softening. The branch softens as a whole for 
q > qAC- The temperature dependences of the sus¬ 
ceptibility at q = 0 (ferroelectric soft mode) and at 
q = 0.1 (TA mode just after avoided-crossing) can 
be described by the Curie-Weiss law with similar 
critical temperatures. 

6 . In Scenario 2 the SETS growth is due to the emerg¬ 
ing central peak. Both TA and TO phonons are un¬ 
altered. The relative rate of growth for the suscep¬ 
tibility, straightforwardly extracted from the cen¬ 
tral peak intensity at q = 0 . 1 , is larger than the 
one for e(T). 

7. Scenarios 1 and 2 are minimal in the sense that 
additional constraints on the fitting parameters do 
not yield plausible fits. The current data does not 
admit a reliable analysis by a more complex model, 
such as Scenario 3. 

2. F-X direction, transverse displacements 

No temperature dependence, no central peak. Excep¬ 
tion - the proximity of the F point. 

3. R-point 

Nearly temperature-independent phonon resonance, 
growing central peak. The relative rate of growth, de¬ 
duced from the central peak temperature dependence, is 
too high to be compatible with the CW law with critical 
temperature below the transition temperature. 

VII. INTERPRETATION AND DISCUSSION 

The aim of the current discussion is mainly limited to 
addressing the question “What kind of phase transition 
can be compatible with the properties of the lattice dy¬ 
namics in the paraelectric phase?”. The non-Bragg scat¬ 
tering experiment in the high-symmetry phase can only 
give the information on the susceptibility of the crystal 
with respect to particular patterns of ionic displacements, 
as defined by Eq. ( 0 ). So we can reformulate the main 


question as “Which trends in the observed temperature 
dependences of different sorts of susceptibility can sug¬ 
gest particular static displacement patterns in the low- 
symmetry phase?”. 


A. Ferroelectricity 

There is no unambiguous evidence whether the corre¬ 
sponding q = 0 soft TO mode should be an underdamped 
phonon or a soft mode of relaxational character. In the 
case of PbZrOs, which can be expected to have dynam¬ 
ics similar to PbHfOs, the infra-red studies show a heav¬ 
ily damped response due to the soft mode^^, supporting 
rather the latter option. Further support comes from 
the results of neutron pair distribution function analysis 
in the high-temperature phase of PbZrOa that indicate 
the presence of quasistatic displacements^^. For PbHfOa 
there are similar, although less detailed evidences^!. The 
study of the dynamics at the zone center was carried out 
by Raman scattering in Ref. [2^ A similar growth of 
the small energy transfer signal for the zone center was 
observed as in the present work for finite wavevectors. 
The authors of Ref. discuss the question of displacive 
versus order-disorder character of the ferroelectric soft 
mode, but do not give an unambiguous conclusion. 


B. Antiferroelectric and incommensurate 
modulations 

1. Displacive-like transition 

The simplest possible interpretation of the IXS spec¬ 
tra in the F-M direction is the following. We admit that 
the model in Scenario 1 is the correct one. This means 
that all the temperature evolution can be explained in 
the soft phonon language. In this case the flat part of the 
TA branch for q > qAc is temperature dependent and the 
rate of softening, daJTA(q) /dT, is similar in magnitude for 
all considered q-points. The data do not reveal any min¬ 
imum of the TA dispersion curve. This is in direct anal¬ 
ogy with the TA branch in PbZrOgSii^. Following Ref. 
one may suggest that the formation of the commensurate 
antiferroelectric modulation occurs due to the first-order 
phase transition conditioned upon Umklapp interactions. 
On the other hand, the incommensurate phase transition 
can also be possible due to the first order transition. The 
absence of a minimum in the dispersion curve, in this con¬ 
text, only shows that the transition should be not very 
similar to the second order phase transition. An adequate 
estimate of this “dissimilarity” could be the difference 
between the critical temperature T™ Ri 295K and the 
actual transition temperature. This difference is about 
170 K. So, the actual transition can be “sufficiently dis¬ 
similar” with the second order phase transition, despite 
the observed mode softening. 
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Z. Order-disorder-like transition 

An alternative interpretation arises when we assume 
that Scenario 2 is correct. The susceptibility Xcp*(q)i 
linked with the central peak intensity via (|5]), increases 
on cooling for the whole T-M direction, so one may expect 
a first-order antiferroelectric or incommensurate transi¬ 
tion with a modulation wavevector along this direction. 
The phonons at q > qac are assumed to be tempera¬ 
ture independent. The quasiharmonic models, such as 
the “polarizability model” from RefJ^, or the rigid ion 
model, as employed by us, suggest that the zone center 
TO phonon shall be also temperature-independent in this 
case. So, the temperature dependences of both the dielec¬ 
tric {q = 0 ) and the generally incommensurate {q > 0 ) 
susceptibilities should be in some way determined by a 
temperature dependence of Xcp*(q). 

First we consider the case in which both these suscepti¬ 
bilities are just proportional to Xcp*(‘1 ) corresponding 
wavevectors. With this assumption it turns out (as was 
already noted in the section dealing with the data treat¬ 
ment according to Scenario 2), that the relative rate of 
the susceptibility increase with respect to modulations 
with q — 0.1 is larger than the one for e(T), at least for 
Tc < T < Tc + AT. On the other hand, Xcp*(^) 
q = 0.1 does not appear to increase slower than the one 
for larger q (as determined in Scenario 2, not shown). So, 
within the present assumption, the relative growth rate 
for the susceptibility at q = 0.1 is larger than the one 
for g = 0 and q > 0.1. This is unlikely. In the presence 
of the ferroelectric soft mode, Xcp*(^) should be a func¬ 
tion, decreasing with q increase^. The assumption that 
the susceptibility is just proportional to Xcp*(q) '''^hh 
respect to modulations is therefore unlikely. 

In an alternative scenario we admit that the TA 
phonon and the central peak are due to displacements 
with similar patterns. This assumption corroborates well 
with the fact that the central peak is observed in the 
same areas of the reciprocal space, where the TA branch 
is unusually flat. In this case it is natural to assume that 
at q = 0.1 and for larger values of q the susceptibility 
with respect to the formation of the modulation waves 
is jointly reflected by the susceptibility with respect to 
the displacements in the TA phonon and by the suscep¬ 
tibility due to the redistribution of the quasistatic dis¬ 
placements. In analogy with that, the dielectric response 
{q = 0) is the result of both the ionic shifts in the TO 
phonon mode and the redistribution of the quasi-static 
displacements (that would be seen as a central peak in a 
scattering experiment). The simplest case assumes that 
the patterns of ionic displacements in the TA branch and 
in the quasistatic disorder are the same. This implies 
equal structure factors for these two components. The 
integral intensity over SETS is thus determined by the 
sum Xcp*'''(q) +XTA‘“(q) = XsETs(q) as (from Eq. (0)) 

^rETs(q) = (13) 


The temperature dependence of ::^ETs(q) i® nearly the 
same as the one, determined for x^*“^(q) in Scenario 1, 
because the parameters of the TO phonon in both sce¬ 
narios are nearly the same implying the same integrals 
over SETS. In this case the consideration of the tempera¬ 
ture dependence of the static susceptibility is completely 
equivalent to the one for the displacive-like scenario. 


3. Generalized viewpoint 

In the preceding paragraphs we considered the two 
viewpoints on the pre-transitional dynamics. In the for¬ 
mer, the crystal approaches the displacive-like transition 
driven by the soft phonon. In the latter, it approaches 
the order-disorder-like transition, where the susceptibil¬ 
ity increase is reflected in the spectra by an increase of 
the displacement-displacement correlation function (cen¬ 
tral peak) for the quasistatic displacements. These two 
viewpoints are the limiting cases of the following more 
general consideration. The dynamic scattering function 
for the lattice mode j is written as^ 

5,(Q,cc) = |G,(Q,q)p^xKq>‘^)> (14) 

where Gj(Q,q) - structure factor for the lattice mode 
and Xj(q, w) - imaginary part of the dynamic suscepti¬ 
bility. In the case of the displacive-like transition view¬ 
point the dynamic susceptibility is the one of the damped 
harmonic oscillator xdho(w, Wj, T) = (w| — — iwT)”^. 

In the case of the order-disorder-like transition, the dy¬ 
namic susceptibility can be represented by xod+dho = 
XDHo(w,a;cp,r®^°) -|-xdho(w, wta, T), where the part 
responsible for the central peak is parametrized by a 
DHO susceptibility function with very large damping 
constant T®^®. The contributions of the two terms to 
the cumulative static susceptibility xIe^ts given by 
the parameters wcp and wta via Eq. ©• Among the 
possible more sophisticated models for the dynamic sus¬ 
ceptibility function we shall mention the one introduced 
in Refs. [2^ and HO and used in Ref. H and d for the inter¬ 
pretation of the IXS spectra in PbZrOa. In the current 
paper we do not attempt using it for PbHfOa because 
the IXS spectra do not contain enough details to make 
its application reliable for practical data analysis. 


4- Summary on the F-M direction 

Prom the considerations above it turns out that it is 
unlikely that we shall associate the central peak in Sce¬ 
nario 2 with the static susceptibility and that the inte¬ 
gral scattering intensity (SETS) should be used for that 
purpose instead. The static susceptibility XsETs(q)’ 
small q, appears to have a similar temperature depen¬ 
dence as the one of e(T). Consequently, the transition to 
the modulated IC phase can be as anticipated from the 
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IXS results, as ferroelectricity can be anticipated from 
the dielectric results. 

Regarding the nature of the dynamic susceptibility, no 
unambiguous conclusion can currently been drawn. Al¬ 
though the analysis of spectra at q = (0.12,0.21,0) sug¬ 
gests that the TA phonon has a temperature dependence, 
the magnitude of this temperature dependence cannot be 
traced quantitatively and the presence of the central peak 
cannot be ruled out. The presence of the central peak is 
natural to expect due to the known evidences in favor 
of the presence of quasistatic displacements in the para- 
electric phase^l. We note that the problem of theoretical 
understanding of the difference between displacive and 
order-disorder regimes of the phase transition, and pos¬ 
sible crossovers between them, is quite challenging, and 
was intensively studied in the context of simple models 
for the Hamiltonian, such as the so-called (j)'^ modeti^. 
The relationship between soft modes and central peaks is 
also an intensively studied topio^. We think that there 
is space for a deeper understanding of the relationship 
between the rather broad spectrum of theoretical results 
and our experimental observations. 


C. R-point order parameter 

At the R-point we observed nearly temperature- 
independent phonons and a strongly temperature depen¬ 
dent central peak. In the absence of conclusive informa¬ 
tion on the displacement patterns in these two spectral 
features we can only sketch the possible variants of the 
interpretation. In the antiferroelectric phase of PbHfOs 
the superstructure reflections at the R-point (in pseu- 
docubic reference frame) are associated with oxygen oc- 
tahedra rotations. We suggest that the growth of the 
central peak on cooling in the paraelectric phase reflects 
the quasi-static disorder in the shifts of oxygen atoms. As 
it was shown in the corresponding data analysis section, 
the temperature dependence of the central peak inten¬ 
sity is unlikely to be consistent with the CW relationship 
with Tq below the transition temperature: the relative 
rate of growth for the corresponding susceptibility is too 
large near the transition. On the other hand, it might be 
possible that this temperature dependence can be con- 
sistentl y d escribed by more sophisticated models such as 
in Ref. [2^. 


D. Comparison with PbZrOs 

The current study reveals both similarities and poten¬ 
tial differences between the characters of the lattice dy¬ 
namics in PbHfOa and PbZrOg^i^. In both crystals there 
is an increase of the static susceptibility with respect to 
the transverse modulations with the wavevector in the P- 
M direction. The measurements in Refs. H andd and in 
the current study were done in different Brillouin zones 
(BZs). In PbZrOa (measurements in the BZ (3 2 0)) 


the central peak was reported for all temperatures up to 
T = 780 K. In PbHfOs (BZ (2 0 0)) there is no cen¬ 
tral peak at high temperatures {T = 773 K) and it can¬ 
not be unambiguously identified at lower temperatures: 
the fits with and without the central peak describe the 
spectra equally well for g-points exactly at the P-M di¬ 
rection. On the other hand, in PbHfOa we recorded the 
well-observed, underdamped up to g ~ 0.2, TO phonon 
branch which apparently has an avoided-crossing with 
the TA branch and is the continuation of the initial trans¬ 
verse acoustic slope at g > gAC- The notable difference 
is in the temperature dependence of the R-point spec¬ 
tra. For PbZrOa the R-point central peak was reported 
to be weakly temperature dependent, while in PbHfOa 
it is strongly temperature dependent near the transition 
temperature. 


VIII. CONCLUSION 

We reported a comprehensive analysis of the inelastic 
X-ray scattering data in PbHfOa. The temperature de¬ 
pendence of the small energy transfer scattering points 
towards the possibility of formation of incommensurate 
or antiferroelectric phases in the course of the first-order 
phase transition. Generally it supports the result of 
Refs. 1^ and i that antiferroelectricity can be a lock-in 
phase in the incommensurate phase transition scenario. 
One might expect the incommensurate phase to occur in 
the temperature range currently assigned to intermediate 
AFE phase(s) in PbHfOa, which corroborates with sug¬ 
gestion of Ref. m, based on the analysis of birefringence 
measurements. 

It turns out that the zone-center ferroelectric soft mode 
(of whatever displacive or order-disorder origin) and the 
generally incommensurate soft mode on the right of the 
avoided-crossing region are likely to correspond to simi¬ 
lar ionic displacements that occur homogeneously for the 
zone-center soft mode and in the form of modulations 
for the generally incommensurate soft mode. This as¬ 
sumption is supported by the fact that the temperature 
dependence of the susceptibility on the left and on the 
right hand sides of the avoided-crossing region is approx¬ 
imately the same. 

The temperature evolution of the IXS spectra in 
PbHfOa, similar to the case of PbZrOa, is quite challeng¬ 
ing for an unambiguous analysis. In this paper we rigor¬ 
ously explored up to which extent one may increase the 
flexibility of the data fitting model without compromis¬ 
ing the conditionality of the data fitting problem. With 
the current energy resolution (3 meV), the problem of 
distinguishing the correct variant of the dynamic suscep¬ 
tibility function appears ill-defined, while the static sus¬ 
ceptibility can be well analyzed. One may expect that an 
improved energy resolution can help settling the question 
of the dynamic susceptibility function in PbHfOa. In this 
context PbHfOa can be more prospective than PbZrOa. 
The spectra obtained in this study for PbHfOa appear. 
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on visual infection, more detailed than the ones from 
Refs. 0 and|a, taken with a better resolution of 1.5 meV. 
In particular the TA phonon peaks are better resolved in 
the current study at high temperatures. 
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Appendix A: Spectra for the F-M direction 

The results for the fits of the data along the T-M di¬ 
rection are shown in Fig. [8] (for q = 0.2 and q = 0.3) and 
in Fig. [9] (for q = 0.4 and q = 0.5). The presentation 
method is the same as in Fig. - for each q there are 
two adjacent panels showing the fits according to the two 
scenarios. Starting from g = 0.3 the signal from the TO 
phonon becomes poorly resolved and appears as an over¬ 
damped DHO signal in the fits. The description of the 
shape is not always satisfactory in this region. In partic¬ 
ular, for q = 0.3 the shape of the experimental curve does 
not strictly follow the fit at energy transfers around 7-8 
meV. The current data quality does not allow a better 
description. 
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Scenario 1, q=0.2 

1=773 K, q=0.2 


Scenario 2, q=0.2 






Scenario 1, q=0.3 

T=773 K, q=0.3 


Scenario 2, q=0.3 

T=773 K, q=0.3 






T=513 K, q=0.3 


T=513 K, q=0.3 




FIG. 8. IXS spectra of PbHfOa corresponding to Q = (2.2, 0.2, 0) and Q = (2.3, 0.3, 0) treated according to Scenarios 1 and 2. 
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Scenario 1, q=0.4 Scenario 2, q=0.4 


Scenario 1, q=0.5 Scenario 2, q=0.5 


T=773 K, q=0.4 


T=773 K, q=0.4 


T=773 K, q=0.5 


T=773 K, q=0.5 












FIG. 9. The IXS spectra corresponding to Q = (2.4, 0.4, 0) and Q = (2.5, 0.5, 0) treated according to the Scenarios 1 and 2. 



















































